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The development of the ducted fan noise propagation and radiation code CDUCT-LaRC at NASA Langley 
Research Center is described. This code calculates the propagation and radiation of given acoustic modes ahead 
of the fan face or aft of the exhaust guide vanes in the inlet or exhaust ducts, respectively. This paper gives a 
description of the modules comprising CDUCT-LaRC. The grid generation module provides automatic creation 
of numerical grids for complex (non-axisymmetric) geometries that include single or multiple pylons. Files for 
performing automatic inviscid mean flow calculations are also generated within this module. The duct propa- 
gation is based on the parabolic approximation theory of R. P. Dougherty. This theory allows the handling of 
complex internal geometries and the ability to study the effect of non-uniform (i.e. circumferentially and axially 
segmented) liners. Finally, the duct radiation module is based on the Ffowcs Williams-Hawkings (FW-H) equa- 
tion with a penetrable data surface. Refraction of sound through the shear layer between the external flow and 
bypass duct flow is included. Results for benchmark annular ducts, as well as other geometries with pylons, are 
presented and compared with available analytical data. 


Introduction 

The prediction of aircraft engine noise, a primary source 
of sound from aircraft, is important for addressing the is- 
sues of community noise and cabin noise control. One 
of the main aims of Computational Aeroacoustics (CAA) 
is the development of physics based methodologies for 
engine noise prediction. This paper describes the devel- 
opment of the ducted fan noise propagation and radiation 
code CDUCT-LaRC at NASA Langley Research Center. 
This code calculates the propagation of a given acoustic 
source ahead of the fan face or aft of the exhaust guide 
vanes in the inlet or exhaust ducts, respectively. Subsequent 
to the propagation calculations, the code has the capability 
of computing the noise radiation field outside the duct. In 
addition to the ducts being acoustically lined in specified 
areas, there are other issues that make the propagation and 
radiation analysis complex: 

1 . The duct may have a non-uniform circular or annular 
section. As a result, the flow in the duct may be non- 
uniform. 

2. There may be a pylon in the exhaust duct resulting in 
one or two C- shaped regions. This will completely 
change the character of acoustic wave propagation 
from that of an annular duct and conventional wave 
propagation theory will then be insufficient. 
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3. The acoustic lining on the duct wall may be non- 
uniform (i.e. circumferentially and radially seg- 
mented). 

CDUCT-LaRC aims to provide a user-friendly environment 
in which to study these and other issues in the pursuit of 
low-noise engine designs. In light of this, the development 
of the code has followed an approach which attempts to 
separate the main functional capabilities into specific mod- 
ules. In this way, as improvements or new approaches be- 
come available, they may be incorporated more efficiently 
and without major modification to the user interface. The 
discussion to follow begins with a description of the ma- 
jor modules within CDUCT-LaRC. This is followed by the 
presentation of various benchmark cases to provide some 
confidence in the capabilities of the code. Preliminary 
results for more complex geometries are then shown to il- 
lustrate some of the areas in which the code may be applied. 
Finally, areas of further development and applicability are 
discussed. 

Modules of CDUCT-LaRC 

CDUCT-LaRC is composed of five distinct modules: in- 
put and output specification, Computational Fluid Dynam- 
ics (CFD) and acoustic grid generation, background flow 
calculation, duct acoustic propagation, and duct acoustic 
radiation. A graphical user interface (GUI) is used for in- 
put and output specification and to guide the user through 
various options offered within the modules. Typically, a 
session begins with the specification or generation of a 
background grid. 

Grid Generation Module 

The most general noise computation, for either an inlet 
or aft-fan case, involves the calculation of a background 
flow (i.e. a CFD calculation). This requires a numerical 
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grid about a possibly complex geometry, the generation of 
which may be outside the interests of many users. In an 
effort to ease this burden, CDUCT-LaRC offers a grid gen- 
eration module. Although a complete description of the 
grid generation techniques is beyond the scope of this pa- 
per, a brief discussion of this module is in order. Following 
guidelines for grid topology, the user simply provides the 
duct surface geometry in PLOT3D surface format. The in- 
put geometry is typically the engine cowl and hub, although 
a by-pass duct might also include one or two pylons. Based 
on the defining topology, a background flow grid is auto- 
matically generated that is constrained to the input surfaces 
and clustered in regions of higher curvature. The user has 
control of grid density, but a complete knowledge of grid 
generation is not required. It should be noted that it is not 
necessary to use the grid generation module to generate a 
background flow grid. It is possible for a user to generate 
and import a CFD grid, in the appropriate PLOT3D format, 
if desired. 

In addition to the background flow grid, an acoustic 
grid is used in the propagation calculations. This second 
grid is necessary because of the conflicting requirements of 
the background flow and acoustic solution methodologies. 
While the background flow grid will necessarily include 
some grid point clustering, this is not desirable for the prop- 
agation calculations. On the contrary, as described subse- 
quently, the propagation calculation benefits from numer- 
ical grids that are as nearly orthogonal and as uniform as 
possible. Within CDUCT-LaRC, this acoustic grid is gen- 
erated directly from the previously created (or imported) 
background flow grid and constrained to the input surfaces. 
Using volume Non-Uniform Rational B-Splines (NURBS) 
techniques, the background grid is splined and the grid 
points are redistributed uniformly to produce a grid that is 
more appropriate for the propagation calculations. As was 
the case with the background flow grid, it is also possible 
for a user to import an acoustic grid if desired. With the re- 
quired grids in hand, the background flow and propagation 
calculations may be performed. 

Background Flow and Propagation Modules 

The background flow module provides the ability to 
study the effects of non-uniform background flow for single 
or multi-block geometries. Examples include engine inlets 
and by-pass ducts containing one or two finite-thickness 
pylons. A well tested and robust flow solver is currently 
used to produce inviscid solutions, although the ability to 
calculate viscous flow fields does exist. The required input 
files are generated within CDUCT-LaRC using supplied 
ambient and duct flow conditions. For situations in which 
it is appropriate, the background flow calculation may be 
avoided and zero or uniform background flow specified. 
As was the case with the background flow grid, a previ- 
ously available flow solution may be imported. Regardless 


of the method in which it is obtained, the background flow 
solution is then interpolated onto the acoustic grid, again 
using volume NURBS techniques. The Mach number dis- 
tribution on the acoustic grid and the input information may 
then be obtained for the propagation calculations. 

The duct propagation module is based on the Boeing 
CDUCT code developed by Dougherty 1 and revised by 
Lan. 2 This utilizes the parabolic approximation and is an 
appropriate initial model in that it can satisfy the three re- 
quirements listed in the introduction. In CDUCT-LaRC, 
the code has been extended to allow a user to automati- 
cally perform multi-block propagation calculations. Within 
this module, the grid connectivity is determined and data 
is transferred from upstream to downstream blocks with- 
out user intervention. Because of the approximations in 
the mathematical model, an orthogonal grid is required for 
the acoustic calculation. However, many geometries will 
require isolated regions in which complete orthogonality 
is not possible. It is believed that such small areas of 
non-orthogonality will not greatly affect solution quality. 
Results of the propagation module include the acoustic po- 
tential or pressure within the duct. If desired, these results 
may be passed to the radiation module for acoustic radia- 
tion calculations. 

Radiation Module 

The duct radiation model is based on the Ffowcs 
Williams-Hawkings (FW-H) equation with a penetrable 
data surface. 3 Testing of this equation for many acoustic 
problems has shown it to generally provide better results 
than the Kirchhoff formula for moving surfaces. Based 
on the background flow conditions and propagation solu- 
tion, this module calculates the radiated acoustic pressure 
at various observer locations. These observer locations are 
specified in the form of a PLOT3D grid file. Currently, the 
data surface is taken to be the inlet or exhaust plane for inlet 
or aft-fan cases, respectively. While this provides excellent 
results in many situations, these choices of data surface lo- 
cation lead to a few limitations. The first appears when 
large aft or forward radiation angles are of interest for inlet 
or aft-fan cases, respectively. In these situations, the effects 
of the engine cowl are not taken into account and the results 
are essentially symmetric about the data surface. The sec- 
ond issue involves the effects of the shear layer on acoustic 
radiation. 

The shear layer between the bypass flow and external 
stream can refract the sound waves radiated to the farfield. 
Radiation results can be improved by including this effect, 
as well as the reflection of the sound in the bypass region 
from the solid surface external to the bypass duct surround- 
ing the core flow (see figure 1). One way to achieve this 
is to extend the bypass duct acoustic calculation, which is 
based on the parabolic approximation, beyond the exhaust 
plane and to apply a boundary condition (BC) similar to a 
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Figure 1 Boundary Regions for the Aft-Fan Geometry 

liner BC on the shear layer. Although a great deal of for- 
malism has been avoided up to this point and the BC is 
not yet implemented, the novelty of this approach makes it 
appropriate to present the boundary condition derivation. 
On the shear layer, satisfy the following two conditions: 


1 . The particle displacement must be continuous across 
the shear layer. 

2. The acoustic pressure must be continuous across the 
shear layer. 


The first result in the derivation involves the acoustic ve- 
locity normal to the shear layer. This follows a derivation 
of the liner BC by Myers 4 as presented by Farassat and 
Dunn. 5 The subscripts 0 and 1 are used for the background 
and acoustic quantities. The particle displacement ampli- 
tude normal to the shear layer is also denoted by q 3 and the 
time dependence is harmonic with the factor e ~ m . Follow- 
ing the notation of Farassat and Dunn 5 [eq. (3) with 8=1], 
one may write 

^ = -mq 3 = [^3no-(«o-Vi^)+i<t-«o] 1 

= [g 3 no-(no-Vi<o)+i<t-«o ] 2 W 


where the superscripts 1 and 2 designate regions 1 and 2 in 
figure 1 , it is the fluid velocity, and it is the unit normal to 
the shear layer. 

From equation 1 , one finds 

<73 = — ~ M _jr — , = — ~ M jr . . (2) 

/co + [no • (no • Vwo)] /to + [no ■ (no ■ Vwo)] 
and 

_ /co + [no-(no-Vwo)] ^ 

u ln ~ . r _). U l rr v 5 ' 

m + [rio-(no-Vuo)] 

This result shows that if u\ n is found on the shear layer from 
the parabolic approximation method in region 1, then u\ n 
on the external side of the shear layer may be calculated. 
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Figure 2 Trace Velocity Description 


Next, the derivation of the impedance, zi, on the side of 
the shear layer adjacent to region 2 is presented beginning 
with the model and assumptions.* First, the shear layer is 
taken to be infinitely thin, an assumption previously uti- 
lized in equation 3. Second, the barber pole pattern of the 
modes present in the flow region is assumed to persist in the 
flow region inside the shear layer. This barber pole pattern 
has a trace velocity on the shear layer surface in the direc- 
tion normal to itself. Figure 2 shows this trace velocity, V^, 
and the method in which it is calculated. This figure also 
shows a local tangent plane to the shear layer surface which 
is cylindrical. Assuming the acoustic pressure in a mode to 
be proportional to exp [— / (co t — m0 — k x x)\ where m is the 
circumferential mode number, v is the axial distance, and 
k x is the axial wave number, it follows that 



m 


x = 


co 

kx 


( 4 ) 


p a °V C 1 p kRs 

R S Q = R s c 1 = ci 

m m 


co/ci k 

±= ^r Cx = k Cy 

Kx K x 


( 5 ) 

( 6 ) 


V s = 



= kc\ 



( 7 ) 


Here, R s is the local radius of the shear layer (assumed al- 
most constant), k = co/ci , and c\ is the local speed of sound 
in region 1 of the shear layer. 

From Morse and and Ingard 6 (equation 11.1.21), the 
impedance, zi, is 


Z2 = 


P2C2 

(1 + Af2 cos (f>) sincf) 


( 8 ) 


*Note that until the end of this section, we work with the zeroth order 
quantities (background flow, time independent quantities). For simplicity 
of notation, we denote regions 1 and 2 by subscripts. 
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Benchmark Results 



\ Plane wavefront 
in region 2 


Figure 3 Effective Wave Speed Description 


where M 2 = 1 * 2/02 is the Mach number based on the back- 
ground flow in region 2 and (]) is the preferred direction 
of propagation of plane waves in region 2. This angle is 
found purely by kinematic considerations. Figure 3 shows 
the geometry in the local plane containing V s and the lo- 
cal normal to the shear layer which is along axis 2. Let 
— (sin cos 4>) be the direction normal to the wavefront. 
The effective wave speed in region 2, ct , is 

ct = C2~^ + W2COS\|/ (9) 


where cos \| f = x/V s is shown in figure 3. The projection of 
c ( t on the shear layer is equal to the trace speed V s . This 
gives 

V s = t/2COS\|/ + C2sin(|) 
from which it is seen that 


sin4> = 


Vs — U2 COS\|/ 


C2 


Using this in equation 8 leads to 

P2C2 (V s ~ U 2 ) 


Z2 : 


V s 


\/(^-«2) 2 


The final step in the formulation is that 


zi 


_ ^ln 
U In 


Zl 


( 10 ) 


(ID 


( 12 ) 


and equation 3 provides an expression for U \J U \n ■ 11 
should again be noted that the radiation results to follow 
use the duct exhaust plane as the data surface and do not 
utilize this boundary condition; however, implementation 
is underway. 
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In order to verify that the various modules function prop- 
erly, a benchmark noise propagation and radiation study 
has been performed using a straight annular duct with mean 
flow. Although the applicability and accuracy of the current 
propagation module has been documented, 2 further confi- 
dence in its multi-block implementation is desirable. Initial 
verification of data transfer and retarded time calculations 
within the radiation module is also of interest. The bench- 
mark geometry is a straight annular duct with dimensions 
that approximate the more realistic bypass duct presented 
in sections to follow. As such, the inner radius is 0.285 m 
(11.2 in), the outer radius is 0.412 m (16.22 in), and the 
overall length is 1.07 m (42.9 in). The mean flow within 
the duct was taken to be uniform in the direction of the 
duct axis at M = 0.4. The specification of uniform mean 
flow meant that a background flow grid and calculation 
were not necessary. Instead, it was possible to simply spec- 
ify the uniform flow conditions on the acoustic grid. As a 
side note, this illustrates some of the flexibility of CDUCT- 
LaRC in that propagation calculations are carried out with 
a previously obtained mean flow solution. 

Propagation Calculations 

The aforementioned geometry and mean flow provide 
an opportunity to perform in-duct propagation calculations 
and formulate analytical solutions for comparison. The 
process begins with results for modal propagation in the 
direction of the mean flow (i.e. positive traveling waves). 
Assuming the duct axis to run in the x-direction from x = 0 
to x = 1 .07 m, this implies modal source specification at 
v = 0 and approximates the situation of an aft-fan case. The 
source frequency was taken to be 5000 Hz and the (10, 1) 
hardwall annular duct mode (corresponding to a cut-off ra- 
tio of 3.59) was specified. Figure 4 shows contours of the 
real part of the complex acoustic potential for this mode 
on a grid of dimension 17 x 113 x 153 (radial, circum- 
ferential, and axial dimensions, respectively). This result 
is compared with an analytic solution obtained from duct 
propagation theory in figures 5 and 6. The plotted data is 
extracted from grid lines at 0 = 0, r = 0.308 m (12.1 in.) 
and r = 0.377 m(14.8in). As these plots show, the solu- 
tions match very well. In order to provide confidence in 
the data transfer between blocks for multi-block propaga- 
tion cases, this problem was also run with the duct split 
into three separate blocks. The result was then compared 
with the single block result presented in figure 4. Figure 7 
provides an example comparison for data extracted along a 
grid line at 0 = 0, r = 0.308 m showing that the single and 
multi-block solutions match. 

In addition to the previous cases, an upstream propaga- 
tion case was also examined. Assuming the duct axis to 
again run in the x-direction from x = 0 to x = 1.07 m, 
upstream propagation implies modal source specification 
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Figure 4 Real Part of the Acoustic Potential - Downstream 
Propagation [M — 0.4, 5000 Hz, (10,1) mode]. 
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Figure 5 Real Part of Acoustic Potential - Downstream Prop- 
agation [r = 0.308 m, 0 = 0,M — 0.4, 5000 Hz, (10, 1) mode]. 


at x = 1 .07 m and approximates the situation of an inlet 
case. The source frequency was again taken to be 5000 Hz 
and the (10, 1) hardwall annular duct mode was specified. 
Figures 8 and 9 present comparisons of the CDUCT-LaRC 
results and the analytic solution from duct propagation the- 
ory. The plotted data is extracted from grid lines at 0 = 0, r 
= 0.308 m (12.1 in.) and r = 0.377 m (14.8 in). As was the 
case with the downstream propagation results, the solutions 
match very well. 


Figure 6 Real Part of Acoustic Potential - Downstream Prop- 
agation [r = 0.377 m, 0 = 0 ,M = 0.4, 5000 Hz, (10, 1) mode]. 



x (m) 


Figure 7 Single and Multi-block Results for the Real Part of 
the Acoustic Potential [r = 0.308 m, 0 = 0 ,M = 0.4, 5000 Hz, 
(10, 1) mode]. 
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Figure 8 Real Part of Acoustic Potential - Upstream Propa- Figure 10 Radiated Sound Pressure Level (dB re: 20 juPa) 
gation [r = 0.308 m, 0 = 0,M = 0.4, 5000 Hz, (10, 1) mode]. 
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Figure 11 Comparison of CDUCT-LaRC and Rayleigh (An- 
alytical) Radiation Sound Pressure Levels [5000 Hz, (10, 1) 
mode]. 


Figure 9 Real Part of Acoustic Potential - Upstream Propa- 
gation [r = 0.377 m, 0 = 0 ,M = 0.4, 5000 Hz, (10, 1) mode]. 

Radiation Calculations 

As mentioned previously, CDUCT-LaRC contains a ra- 
diation module based on the Ffowcs-Williams Hawkings 
equation with a penetrable data surface. The benchmark 
study also entailed using the benchmark aft-fan (down- 
stream) propagation results to calculate the subsequent 
noise radiation field outside the duct. Radiation calcula- 
tions were performed on a hemisphere of radius 1.91 m 
(75 in) centered on the duct axis in the exhaust plane. Fig- 
ure 10 shows the radiated pressure (dB re: 20 /iPa) on this 
hemisphere. Since the benchmark geometry is a straight 


annular duct, a Rayleigh piston in a wall formulation may 
be used to obtain an analytical result for the exhaust ra- 
diation. Figure 1 1 shows the comparison of the CDUCT- 
LaRC and Rayleigh results along an arc on the observer 
hemisphere. The solutions compare well as the relative 
amplitudes and directivity pattern are similar with slight 
discrepancies in the overall amplitude of some of the mi- 
nor lobes. In addition to the (10, 1) mode, the full series of 
propagation and radiation calculations were carried out for 
the (10, 4) mode. For the benchmark duct and source fre- 
quency of 5000 Hz, this provides a situation in which the 
propagating wave is much closer to cut-off (cut-off ratio of 
1.26). Comparison of the CDUCT-LaRC and Rayleigh re- 
sults along an arc on the observer hemisphere is presented 
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Figure 12 Comparison of CDUCT-LaRC and Rayleigh (An- 
alytical) Radiation Sound Pressure Levels [5000 Hz, (10, 4) 
mode]. 

in figure 12. The relative amplitudes and directivity pattern 
again compare well and the results of these various bench- 
mark problems provide some justification for performing 
preliminary calculations on more complex geometries. 

Preliminary Results 

In an effort to incrementally work toward more complex 
geometries, some preliminary calculations were performed 
on two ducts similar in dimension to the benchmark straight 
annular duct. The first of these ducts (duct 1) is identical 
to the benchmark duct except that an infinitely thin pylon 
is located in the middle third of the duct. The second is a 
by-pass duct (duct 2) similar to that of a business jet. Both 
ducts have the same inlet geometry and overall length of 
1.07 m (42.9 in), however, duct 2 has non-constant radial 
dimensions. Additionally, both ducts have pylons in the 
same axial location with the pylon of duct 2 being mod- 
eled as a NACA 0015 airfoil. Prior to the propagation 
and radiation calculations, the background flow within the 
ducts must be established. Since duct 1 is a straight an- 
nular duct and inviscid flow is assumed, the background 
flow was taken to be uniform at Mach number 0.4 (as with 
the benchmark cases). The background flow module of 
CDUCT-LaRC was used to obtain the mean flow within 
duct 2 with an inlet Mach number of 0.4. 

With the background flow specified, propagation calcu- 
lations were performed within the two ducts. Acoustic 
source specification followed that of the benchmark cases 
and consisted of specifying the (10, 1) hardwall annular 
duct mode at the duct inlet with a source frequency of 
5000 Hz. Figures 13 and 14 show contours of the real part 
of the acoustic potential for ducts 1 and 2, respectively. In 
these figures, the fluid flow and modal propagation are in 
the positive v-direction (i.e. the duct exhaust plane is in the 
foreground). While the typical modal pattern is present, the 



Figure 13 Real Part of the Acoustic Potential - Duct 1 
[5000 Hz, (10,1) mode]. 



Figure 14 Real Part of the Acoustic Potential - Duct 2 
[5000 Hz, (10,1) mode]. 

effects of the pylon are visible as standing waves are gen- 
erated within the pylon region. The differences are more 
evident when the cases are compared with the benchmark 
annular duct contours shown in figure 4. 

The change in modal content and non-uniform mean 
flow resulting from the pylons may also modify the radi- 
ated sound field. As such, the propagation results for ducts 
1 and 2 were used to perform farfield radiation calcula- 
tions. Using the observer grid of the benchmark cases, 
radiation results were obtained on a hemisphere of radius 
1.91 m (75 in) centered on the duct axis in the exhaust 
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Figure 15 Radiated Sound Pressure Level - Duct 1 (dB re: 
20 ^Pa) 



Figure 16 Radiated Sound Pressure Level - Duct 2 (dB re: 
20 A/Pa) 

plane. The radiated pressure (dB re: 20 A/Pa) on this hemi- 
sphere for ducts 1 and 2 are presented in figures 15 and 16, 
respectively. Comparison of these radiated pressure pat- 
terns with the benchmark case (see figure 11) clearly show 
a modification in the directivity pattern. In light of these 
results, a focused study of the effects of pylon placement 
and geometry may prove to be a fruitful area of application 
for CDUCT-LaRC. Additionally, the ability to model more 
complex, multi-block geometries may widen the scope of 
inlet or aft-fan noise calculations that may be investigated. 


Concluding Remarks 

The development of CDUCT-LaRC has followed a mod- 
ular approach with a main goal of providing a physics 
based acoustic propagation and radiation tool within a 
user-friendly environment. Compartmentalizing the main 
functionality provides the ability to make modifications 
or replace entire modules efficiently and without major 
changes to the user interface. The current capabilities of the 
grid generation, background flow calculation, duct acoustic 
propagation, and acoustic radiation modules which make 
up CDUCT-LaRC have been described. With the imple- 
mentation of the shear layer BC formulated herein, larger 
propagation angles will be attainable and the refraction ef- 
fects of the shear layer more accurately modeled. It should 
also be noted that while the benchmark and preliminary 
results focused on hardwall ducts, the ability to handle ra- 
dially and circumferentially varying liners exists within the 
propagation module. 

The presented results provide some indication of the ca- 
pabilities of CDUCT-LaRC, as well as research areas in 
which the code may be applicable. In particular, the modi- 
fication of the modal content by the presence of the pylons 
in the duct, as seen in the results for ducts 1 and 2, is of par- 
ticular interest. In light of this, implementation of modal 
decomposition in the post-processing functionality of the 
code has been initiated and is near completion. Addition- 
ally, the ability to handle complex, multi-block geometries 
within CDUCT-LaRC allows investigation into the effects 
of design modifications such as radial or circumferential 
splitters. While there is a great deal of testing to be per- 
formed, preliminary results show that the modules for duct 
propagation and radiation offer appropriate initial models 
for noise prediction. 
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